Dynamics of inertial vortices in multi-component Bose-Einstein condensates 
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With use of the nonlinear Schrodinger (or Gross-Pitaevskii) equation with strong repulsive cubic 
nonlinearity, dynamics of multi-component Bose-Einstein condensates (BECs) with a harmonic trap 
in 2 dimensions is investigated beyond the Thomas- Fermi regime. In the case when each component 
has a single vortex, we obtain an effective nonlinear dynamics for vortex cores (particles). The 
particles here acquire the inertia, in marked contrast to the standard theory of point vortices widely 
known in the usual hydrodynamics. The effective dynamics is equivalent to that of charged particles 
under a strong spring force and in the presence of Lorentz force with the uniform magnetic field. The 
inter-particle (vortex- vortex) interaction is singularly-repulsive and short-ranged with its magnitude 
decreasing with increasing distance of the center of mass from the trapping center. "Chaos in 
the three-body problem" in the three vortices system can be seen, which is not expected in the 
corresponding point vortices without inertia in 2 dimensions. 

PACS numbers: 03.75.-b, 05.45.-a, 05.60.Gg. 



I. INTRODUCTION 

Recently, there has been much interest in theoreti- 
cal and experimental studies on trapped atomic Bose- 
Einstein condensates (BECs) [Jll] . The superfluid prop- 
erty of atomic BECs arises from a dual aspect of waves 
and particles, i.e., matter waves, and is theoretically de- 
scribed by the macroscopic wave function. Because of the 
nonlinearity of the system caused by interaction between 
particles, the macroscopic wave function can take a form 
of various solitons such as bright, dark, grey, and vortex 
solitons, and these solitons are experimentally observed 
in BECs. 

As for bright solitons, Martin et al.^ theoretically pre- 
dicted that the particle-like behavior of three bright soli- 
tons in a one-dimensional 87 Rb BEC was non-integrable 
and showed its change from regular motions to chaos. 
Perez-Garci'a et ah Q applied a variational method to dy- 
namics of bright solitons in a two-dimensional (2-d) BEC 
and showed that the center of mass of each soliton obeys 
Newtonian dynamics and Ehrenfest's theorem is valid if 
the phase of BEC wave function will be suitably chosen. 

As more interesting systems, we can consider soli- 
tons in multi-component BECs which consist of different 
kinds of atoms or same kinds of atoms having different 
spin and have been experimentally realized. In multi- 
component BECs, there are not only intra-component 
particle interaction but also inter-component particle in- 



teraction which is another origin of nonlinearity, so we 
expect novel soliton dynamics which is not seen in single- 
component BECs. Yamasaki et al.\(^ developed a varia- 
tional method to describe bright soliton dynamics in 2-d 
multi-components BECs, and proposed a model of con- 
servative chaos. In 2-d and 3-d systems, however, bright 
soliton are unstable unless intra-component interaction 
oscillates between attraction and repulsion or coexist- 
ing intra-component quintic (three-body) interaction is 
strong enough. 

On the other hand, topological vortices known as quan- 
tized vortices, i.e., topological defects of the macroscopic 
wave function can be stable in two dimensions. Vor- 
tices in single-component and multi-component 
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12] BECs have been realized experimentally, giving a 



good candidate to study the dynamics of vortices in 2-d 
and 3-d BECs. Most of the theoretical studies, how- 
ever, are limited to the Thomas-Fermi regime (TFR) in 
a single-component BEC 0, • While dynamics of 

the macroscopic wave function of BEC is described by 
Gross-Pitaevskii equation (GPE) in Eq.© below, TFR 
suppresses the kinetic energy part in GPE in the lowest 
approximation and therefore the lowest-order wave func- 
tion cannot have a healing length which is a hallmark of 
the vortex core. It is highly desirable to construct the 
effective theory of vortices beyond the TFR. 

In this paper, we consider the vortices in 2-d multi- 
component BECs beyond the Thomas-Fermi regime. To 
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consider the effective dynamics of point-like vortices, we 
extract some degrees of freedom by using the variational 
approach, and derive an effective dynamics with finite 
degrees of freedom. For the case of a single-component 
BEC without a trap, it is well known that the effective 
hamiltonian for many vortex systems has a standard form 



(1) 
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in the limit of infinitesimal vortex cores. Here m is 
the winding number of the z-th vortex, and = 
y/ (xi — Xj) 2 + (i/i — yj) 2 is the distance between cores of 
i-th and j-th vortices. Equation (TJJ shows that there is 
no momentum degree of freedom, and coordinates Xi and 
yi are formally conjugate each other. On the other hand, 
in the case of multi-component BECs in a trap with each 
component having a single vortex, we shall see vortices 
acquire momentum degrees of freedom or inertia. A very 
preliminary idea of the pre sent work was reported in the 
conference proceedings [l9j. 

This paper is organized as follows. In Sec. UTJ starting 
from the multi-component GPE, we apply the variational 
method with use of vortex solitons in the Pade approxi- 
mation, and derive the effective Hamiltonian for vortices. 
Confining to the case of two and three vortices, we nu- 
merically investigate the detailed dynamics of vortices in 
Sec. IIII1 There we see a chaotic behavior of the system 
consisting of three vortices. Scction llVl is devoted to con- 
clusions and discussions. 



II. EFFECTIVE NONLINEAR DYNAMICS 
GENERATED BY THE MULTI-COMPONENT 
GPE 

In this Section, we consider the trapped multi- 
component GPE with vortices and extract some degrees 
of freedom of vortex soliton by using a variational tech- 
nique. 

BEC at zero temperature is described by 
the GPE. We shall consider a 2-d system of 
trapped n-component macroscopic wave function 
$ 1 (t,x,y),$ 2 (.t,x,y),- ■ ■ ,$ n (t,x,y) satisfying the 
equations 
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-V 2 + V(x,y) + gii \^i{t,x,y)\' 



*i(t,x,y), (2) 



for i,j — l,...,n. Here the normalization condition 
for each component of wave functions is defined by 
J \$i(t, x, y)\ 2 dxdy = 1 after a proper rescaling of $j by 
the particle number N common to all components. The 
effect of trapping is expressed by V(x,y) = (x 2 + y 2 ). 
Equation ^ is expressed with use of scaled variables: 
using the confining length I = \J and oscillation pe- 
riod t = ui , space coordinates are scaled by I, time by 
2t, wave function by j, and nonlinearity by j— . The 
nonlinearity coefficients, gty = SffNaij/l with ay the 
scattering length for binary collisions, are assumed to 
be positive and much larger than unity, gu and gij with 
i 7^ j stand for intra-component and inter-component in- 
teractions, respectively. We shall choose gu = <7i(^> 1) 
for all i and gij = 32 1) for i ^ j. 

In the absence of the inter-component interaction, each 
component has stationary states of a vortex. So, we con- 
sider the case in which each component has one vortex 
and vortices interact with each other through the inter- 
component interaction. Our goal is to derive from (j2Jl the 
evolution equation for the collective coordinates of trial 
vortex functions (TVFs). The collective coordinates for 
a vortex are phase variables besides the coordinates of 
a vortex core. We shall use TVF beyond the Thomas- 
Fermi regime, by incorporating the effect of a kinetic en- 
ergy in GPE in Eq.©: As for the amplitude of TVF, 
we choose a vortex function based on the Pade approx- 
imation [2^, EH which is regularized due to a trap. As 
for its phase, we Taylor-expand the phase with respect 
to space coordinates around the vortex core. Then TVF 
with winding number rij = ±1 is given by 



$i(t,x,y) = fi(t,x,y)exp[i<j>i(t,x,y)] 
1 -1- r 2 

N exp 



x z + y l 



2A 



(x - Xi) 2 + (y - y t ) 2 



2£ 2 + (x - Xi ) 2 + {y - yi f 



x exp 



n, tan 



1 V-Vi 



+ a>i(x - x^ + /3i(y - y t ) 



(3) 



with the normalization factor N = 

00 

=4 ■ Here T(0,z) = J t^e^dt 

is the incomplete gamma function of the second kind, 
whose expansion with respect to z is given in Eq. (|A10p 
in Appendix [A"l 

The collective coordinates are locations of the core 
(xi,yi) and the first-order coefficients (a^ft) of Taylor- 
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expansion of the phase 4>i(t, x,y) with respect to (x — 
x iiD — Di)- A in the Gaussian amplitude factor re- 
flects a trap in Eq.©. £ is the healing length re- 
lated to vortex core size. The condition to minimize 
the energy E = J dxdy (|V$,| 2 + V(x, y)\^\ 2 + f |^| 4 ) 
for the individual static component in Eq.© centered 
at the origin, leads to A 

KJttV* -1/4 



fi- and 



respectively, where j(= 0.57721) is 



Euler constant. We shall use these values for A and £ in 
this paper. 

The form in Eq.©, which is a product of the vortex 
solution in the absence of a harmonic trap and Gaussian 
factor due to the trap, gives a suitable TVF for a vor- 
tex under the strong nonlinearity. A different form with 
use of eigenstates (with non-zero angular momenta) un- 



der the 2-d harmonic trap 22H24| has no small healing 



length, results in the inter- vortices force growing with 
inter-vortices distance, etc, and cannot be suitable as 
TVF under the strong nonlinearity. 

First of all we note: GPE in Eq.© can be derived 
from the variational principle that minimizes the action 
obtained from Lagrangian density C for field variables, 



£ =E 



> 2 + y 2 )l<J> l | 2 + |l^l 4 ]+E^I^I 2 l^l 



In fact, Eq.© is obtained from Lagrange equation: 
d dC 



dt d$* 9$ 



dC dC 



<9V$* 



0. 



(4) 



(5) 



We now insert TVF in Eq. © into Eq. ©. Noting 
{xi,yi) and (ai,Pi) as time-dependent variables, Eq. © 
becomes 



C = 





+ f/f+^+y 2 )/ 2 



By integrating £ over space coordinates (x, y), we ob- 
tain the effective Lagrangian L for the collective coordi- 
nates: 



L = 



dxdyC. 



(7) 



In the limit of £ 2 /A <C 1, L is expressed by 



L 



-) E 

N I Z - 



V (xiUi - ViXi) e * I di + -^- L 



+ (a 2 + /3 2 - {oiiXi + PiVi)) d 
- (diXi + fayA 



I 2 



i _ii 

A 



e~* [ 2n 2 d 3 + d A + 2n 2 4 5 



" l A 



n- i 2 / I 2 

+ 2— {yiOi - XiPi) e"* I d 1 + -±-d 2 



+ e & (d 8 + l 2 d 9 ) 



No) ttA U +are 



E n 'v;n- 



(8) 
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Here, k — \J x 2 + y 2 and N — 7/=^- Expressions (: 
Taylor-expansions with respect to £ 2 / 'A) for coefficients 
d\ ~ (ig are listed in Table 1 and the derivation of typical 
coefficients is described in Appendix [A] The expression 
for the inter- vortices interaction U(rij,l % o) , which is a 



function of inter- vortices distance Tu = 



distance of the center of masses from the origin Iq = 



(ig) 2 + (yjj) 2 ' wm ^ e gi ven m Appendix [Bl 
Lagrange equations of motion for the phase variables 
Ui and Pi, 



-( 


f dL\ 


dL 






doti 




(<*) 


dL 


- 

dt 




~dp~ 







(9) 



lead to 



a.i ~ Siajj - rnB 2 yi 
Pi ~ Bij/j + n i B 2 x l . 



(10) 



Equation (flU)) shows that (a,, ft) correspond to gener- 
alized momentum conjugate to (xi,yi) under the vector 
potential. Then the equation of motion for [xi,y{) 




dL 

dxi 

dL 







0, 



(11) 
(12) 



combined with Eq. flTUl) . gives 

x'i ~ niB 3 yi - B 4 Xi - B 5 du ^; VijS> 

I 

Table 1. Expressions for coefficients in Eq.© (: Taylor-expansions with respect to £ 2 /A) . 



Coefficients 


Expressions 


di 


1 + Ml + Mi ln Ml 

1 T A A A 


d 2 


I _ £ 


d 3 


2 2 111 A ' A A A A 




1 + li + Ml + ^iinMl 

2 ^ A ^ A 'A A 


d 5 


1 + 2«1 + £ In Ml 

2 T A ^ A A 


d 6 


1 + Ml + 62£l + Mli n Ml 

1 ' A ^ A ' A A 


d 7 


A A ' A A 


d$ 


A-2£ 2 


dg 


1+Ml + Ml + MlhiMl 

1 ^ A ' A ' A A 



Coefficients B\ ~ B5 are given by 



di - 1 + ^ (7 + l)? 2 



Si 



2di 



Bo = — 



B, 



\d\ 



Mi 

A 



B 4 = 



B 5 = 



A(di - 1 
4didg 

(d 1 -l+ 2 #) 2 " (7 + l) 2 C 4 



(7 + 1) 2 ? 4 



A 2 



N 



2di 



(di 



Mi 

A 



2(7 + l) 2 C 4 



(14) 



Equations (fT0|) and (|T3|) are valid aside from a mul- 
tiplicative global factor (1 + 0(4-)). The smallness of 



will be justified a posteriori. Equation (|13j) shows 
that dynamics of coordinates (xi,yi) is very similar to 
charged particles with charges rij = ±1 under a strong 
spring force with a force constant B4 — 0(<? 2 ) and in 
the presence of Lorentz force with the magnetic field 
B = (0, 0, B 3 ). One should note that the spring force here 
has nothing to do with that of the original harmonic po- 
tential V(x, y) in Eq.Q. The Hamiltonian corresponding 
to Eqs. (fT3j) can be given by 



H 



1 



^ L2m 



[ Pl - n.Ai} 2 + Wi 



(15) 



with use of the momentum p i = mr + rijAj, unit mass 
to = 1, the vector potential Aj = ^-(—yi,Xi), and the 
scalar potential Wi — i-B 4 (x 2 + yf). 

All values Bi ~ i? 5 depend on the strength of inter- 
action gi which is tunable by Feshbach resonance. In 
the effective particle dynamics described by Eqs. (fT5)) and 
(p~5|) , the spring constant B4 is large enough to guarantee 
particles with unit mass to be confined in the neighbor- 
hood of the origin, i.e., the trapping center. This finding 
i 2 

justifies ^ < 1 and will also be utilized to obtain the 
inter- vortices interaction in Appendix FBI 

The scaled inter- vortices interaction U(rij, Iq) is given 
with use of Eq. (|B10[) as 



C7(r y ,/g) = B b U 4 (r l3 ,l%) 




4(^) 2 ' 



2A 



with the coupling constant 
G = 



(16) 



(17) 



( 7 +l) 2 7T 1 /2- 

Therefore the inter-vortices interaction is singularly- 
repulsive and short-ranged with respect to ry with its 
magnitude decreasing with increasing Iq . 

Compared to Eq. (|Tj) , it is clear that the system has 
momentum degrees of freedom, and vortices have a be- 
havior of particle with inertia rather than that of vortex 
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point without inertia widel y us ed in the conventional the- 
ory of hydrodynamics 25l-l29j|. This is one of the main 
assertions of the present paper. The inertia of vortex ap- 
pears already in the single-component BEC with a trap. 

In closing this Section we should comment: we also at- 
tempted to apply the collective coordinate method using 
a Laguerre-type trial function which is a good candidate 
for TVF in the case of a weak nonlinearity 22h24| . How- 



ever, such TVF in the case of a strong nonlinearity proved 
to result in: 1) a time-dependent mass for each vortex 
and 2) the inter- vortex interaction growing with increas- 
ing the inter-vortex distance. Hence this TVF was not 
suitable to describe a dynamics of vortices in BECs with 
a strong nonlinearity. On the other hand, the effective 
vortex dynamics in the Thomas-Fermi regime 
which suppresses the kinetic energy in constructing TVF, 
leads to neither non-zero inertia nor Lorcntz force. 



III. DYNAMICS OF TWO AND THREE 
VORTICES WITH INERTIA 




We shall now focus on the system of two vortices with 
equal winding numbers, and see how the trajectory gen- 
erated by effective particle dynamics in Eq. (fT3)) and ([15)) 
well mimics the orbit of the singular points of wave vor- 
tices calculated by using GPE in Eq.Q. We shall then 
move to the system of three vortices with equal wind- 
ing numbers, and find that chaos appears even in the 
three vortex system. This feature is different from that 
of point vortices system in a single component BEC in 
which chaos can appear in the case of more than three 
vortices. 



FIG. 1: Two vortices dynamics with identical winding num- 
bers (ni = ri2 = 1). gi — <?2 = 100. Initial values: 
xi (0) = yi(0) = -x a (0) = -j/ 2 (0) = ^; xi(0) = aja(0) = 0; 
y'l (0) = —2/2(0) = 1.1. Solid line is a trajectory for relative 
coordinates of a pair of particles calculated by iterating the 
canonical equations for the Hamiltonian (|18[) : crosses are cor- 
responding results for a pair of singular points of interacting 
vortices constructed from the numerical iteration of GPE in 
Eq.©; four sub-panels are wave function patterns for inter- 
acting vortices. 



A. Dynamics of two vortices 

For two vortices with the same winding numbers m = 
712 = 1, Hamiltonian (|15|) can be rewritten as 

H = \(pIt+P 2 v ,t+pIr+pI,r) 

~ "Y (Px,TVT - Py,TX T + Px,RVR ~ P v ,RXr) 

+ + (4 + 4 + 4 + 4) 

- 2 G 2 , exp f_ 4 + 4 + 4+4 ^) i(lg) 



2(4 + 2/1?.) 



A 



where (xt,Ut) and { x r,Ur) play the role of the 
center-of-mass and relative coordinates, respectively and 



(Px,t,Pv,t) and (p x ,R,Py,R.) are their canonical-conjugate 
variables. To be explicit, 

{x t ,Vt) = ^(£1+^2,2/1 + 2/2) 

(xr,vr) = -j=(x 1 -X2,y 1 -y 2 ) 

(Px,t,P v ,t) = ■^=(Px,l+Px,2,Pv,l+Py,2) 

{Px,R,Py,R) = -^(Px,l-Px,2,Py,l-Py,2)- (19) 

Hamiltonian (|15|) cannot be reduced to two inde- 
pendent two degree-of- freedom subsystems (: center-of- 
mass system and relative-coordinates system) because 
the inter- vortex interaction depends not only on but 
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also on Iq. However, we see mostly KAM tori in this 
system. As shown in Fig. [TJ we find that the trajectory 
of {xr,ub) generated by effective two-particle dynamics 
well mimics the corresponding orbit obtained by a pair of 
singular points calculated by using GPE in Eq.(f2]). This 
fact justifies the validity of our trial function in Eq.([3j 
and the resultant equation of motion for collective coor- 
dinates in Ea.(|T3l) and (|T5]). 

B. Dynamics of three vortices 

Encouraged by the effectiveness of the collective coor- 
dinate method in the case of two vortices, we proceed to 
the dynamics of three vortices with the identical wind- 
ing numbers 711=712 = ^3 = 1, whose Hamiltonian (|15p 
becomes 

Tj 1/2 i2 ,2 i2 ,2 ,2\ 

H = ^^C+PyC+PxT+PyT+PxR+PyR) 



- -^{ x CPyC + X R p yR + X T PyT 
~ VCPxC ~ VRPxR ~ VTPxT) 




where Jacobi coordinates (xT,VT,etc) are defined by 



(xt,Vt) 


= Tz ixi 4 


- X 2 + ^3,2/1 + V2 - 




(xcyc) 


= 7=2 {X1 " 


- x 3 ,yi - 2/3), 




(x R ,Vr) 


= T6 ixi 4 


-x 3 - 2x 2 , yi + y 3 


-2y 2 ) 


(PxR,P v r) 




+ Px3 - 2p x2 ,Pyl - 


Vp V 3 - 2p y2 ) 


(PxT,PyT) 


= T3 iPxl 


J rPx2 + Px3,Pyl + 


Py2+Py3) 


(PxCPyc) 


= T2 iPx2 


-Pxl,Py2 -Pyl), 


(21) 



which conserve the canonical structure ({xt,Pxt} — 
3 

§ Y,{ x hPxj) = !> etc -)- Here (xT,yr): (xc,Vc), and 

(xr, yR) represent the center of mass of all three com- 
ponents, the relative displacement, and the bisector of 
the vertex (2:2,2/2), respectively. In marked contrast to 
the massless three vortex system in 2 dimensions which 



is integrable, all 6 degrees of freedom are coupled and 
the number of independent constants of motion is 2 (en- 
ergy and z-component of angular momentum). Then 
Poincare-Bendixon's theorem guaranties the nonintegra- 
bility and chaos of the three inertial vortex system [30j ■ 

We can construct from (|20|1 the canonical equa- 
tions of motion for XT,yT,xc,yc, x R,yR and their 
canonical-conjugate variables, which is solved numeri- 
cally. Poincare cross section and power spectra for the 
kinetic energy Ek (t) in Fig. [2] give a clear evidence of 
high-dimensional chaos. Because of the short-range na- 
ture of the interaction, we see the emergence of chaos in 
low-lying energy regions where vortices often meet each 
other. 




-1 \ 

0.1 0.2 0.3 0.4 
CO 



FIG. 2: Three vortices dynamics with identical winding num- 
bers (n,i = TI2 = W3 = 1). gi = g-2 — 100. Initial values: 
x R (0) = xc(0) = x T (0) = 1, y R (0) = -1, yc(0) = yr(0) = 1; 
Px,r(0) = p x ,c(0) = Px,t(0) = p y , R (0) = Pv,c(0) = Pv,t(0) = 
1. (a) Poincare cross section; (b) Power spectra for kinetic 
energy E k (t). 
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IV. CONCLUSION 

We explored vortex dynamics in the 2-dimensional 
multi-component BEC in the harmonic trap in the case 
that each component has a single vortex. The inves- 
tigation beyond the Thomas-Fermi regime is made on 
the nonlinear Schrodinger equation with strong repul- 
sive cubic nonlinearity. With use of a trial vortex func- 
tion based on the Pade approximation which is regular- 
ized due to a trap, we applied a collective coordinates 
method, obtaining an effective nonlinear dynamics for 
vortex cores (particles), which is equivalent to charged 
particles with inertia under a strong spring force and in 
the presence of Lorentz force with the uniform magnetic 
field. The inter-particle interaction is singularly-repulsive 
and short-ranged with its magnitude decreasing with in- 
creasing distance of the center of mass from the trapping 
center. The most important finding is the nonzero iner- 
tia of vortices, which is not present in the conventional 
theory of point vortices widely used in the usual hydrody- 
namics 25-2t|. The system of three vortices with inertia 



can be chaotic, in contrast to the corresponding case of 
point vortices without inertia. 
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Appendix A: Calculation of typical integrals 



We shall calculate some integrals used in this paper. 



1. Integral A = J J dxdyfi + 
Substituting Eq.Q into the above integrand, we find 



oo oo 



.4 



^ (U i i dxd y^ 

— oo — oo 

xiiy - yi) - y\{x - xj) 

i 2^+( X _ Xl )2 + (j / „y 4 )2 

(x - Xi) 2 + (y - Vl ) 2 



(atiXi + Piiii) 



2£ 2 + {x- Xif + (y - Vl ) 2 



(Al) 



where Nn = , . and ./V is the normalization constant 
defined below Eq.Q. Equation (|Al[) is a sum of n r 
dependent term (Ai) and (ct^, /3i)-dependent term (^2)- 



Below we shall concentrate on Ai because A2 is easily 
calculable. Using polar coordinates as 



x — Xi — p cos 9 Xi = U cos fa 
V - Vi = P sin 6» y t = U sin fa 
dxdy = pdpdO 

we rewrite the integral A\ as 



rii ( N\ t p 2 +<? p 2 



(A2) 



2tt 



x / dO(xi sin 6* — y\ cos#)e~ 



■2plj 



(A3) 



The ^-integration gives 



Z7T 

J d9(xi sin 6 — y\ cos 9)e~ 



2£i cos 



d9 



(x cos fa + y sin fa ) sin ( 



+ (x sin fa — y cos fa) cos 9 



2pl i 



cos 



= -n(xiyi - y t x t ) [ — + 



(A4) 



Here we took 9 = 9 — fa and used the formulas like 

-z cos 9 



d9 cos 9e 



-2tt/i(«), 



(A5) 







where Ii(z) is the modified Bessel function, which is ex- 
panded as 

~3 



As for p-intcgration, we use identities like 

2£ 2 + p 2 9 2i 2 + p 2 
and then apply the formulas 



, r, " 2 1 A 71 + 1 + 1 x 

dpp n e-^ = -A~T( — - — ) 



(A6) 



(A7) 



(A8) 



and 




2e+ P 2 2 e r A 



(A9) 
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r(0,z) = ft 1 e t dt above is the incomplete gamma 

z 

function of the second kind and can be expanded as 
r(0,z) = - 7 -mz + 2-^ + ^ + O(z 4 ). (A10) 
The final result is 

a (N\ 2 n if . . ^_d 2 l 2 \ 
ViVo J A ' ~ VlXl ' e V ~~A~ ) ' 

(All) 

where d\ and d,2 are listed in Table 1. Equation (|A11I) . 
combined with 

A 2 = -(J^) (oiiXi + fcy^d!, (A12) 

gives rise to the calculated result for A in Eq. (|Aip . 



2. Integral B=^JJ ffdxdy 



D 



91 ( N 



2tt 2 A 2 V A 



2(x^+y' ! ) 

dxdye A 



(x - Xi) 2 + (y- in) 2 



M 2 + (x - Xl ) 2 + (y - Vl ) 2 l ' (VI>!1 

Using polar coordinates in (IA2|) . we rewrite the integral 
as 



B = [ !»!?>< ' 



2tt 2 A 2 V A 



o 
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x I _.f - \ I dih (A 14) 



2£ 2 + p 2 



o 



First, we carry out the 9- integration. Note that 



2tt I {z) 



(A15) 



where Io(z) is the modified Bessel function, which is ex- 
panded as 

J (z) = l + ^ + ^+O(z 6 ). (A16) 

Concerning the p-integration, we first employ the de- 
composition 

p 4 4f 2 4£ 4 

1 - ... „ + (A17) 



(2C 2 + P : 



2^ + P 2 (2^ + p 2 ) 2 ' 



and then apply the formulas like 



_ 2 P - i Apk \ A 2'? 
dppe -=-J [ — 1 = jeT. 



(A18) 



The final result is 



b = i^V4 

V A y tta 



- + d 7 e~^ 
4 



(A19) 



where d 7 is given in Table 1. 

Appendix B: Calculation of inter- vortices 
interaction U 

This interaction is due to the integral 

oo oo 

U = 92 J J tfffdxdy 

— oo — oo 

^ oo oo 

^ i . / dxdye~ 



— oo — oo 
2 i /„, „, \2 



7T 2 A 2 ^^0 

(a - Xjf + (y- yi) 
2£ 2 + (x - Xi) 2 + (y- Vi ) 2 

{x - x 3 ) 2 + {y - yj ) 2 
2e + (x-x 3 ) 2 + (y- yj ) 2 - 



(Bl) 



With use of a prescription U = ^i 2 A ? (j^j U, the inte- 
gral becomes: 



U 



dxdye 



1 - 2£ 2 



1 



2£ 2 + (x-x,) 2 + {y- yi ) 2 
4£ 4 



+ (*->■ 3) 



(2e + (x-x l ) 2 + (y-y l ) 2 )( i^j ) 
*7 + U 2 + U A . 



(B2) 



Uq, U2 and K4 corresponds to the contributions from 
0{£°), 0(£ 2 ) and 0(£ 4 ), respectively Among them, U 4 
is responsible to the vortex- vortex interaction, which we 
shall calculate below. 

Let's define the center-of-mass and relative coordinates 
by 



Xj + Xi ij ^ yj + yi 



(B3) 



9 



2tt 




FIG. 3: New integration variables, center-of-mass and relative 
coordinates. 



and 



Xi Xj — Tij COS (pij , JJi yj — V^j SI 



(B4) 



respectively, and transform the integration variables to 
new ones as (see FigEJ) 

x — x 1 q — pcos9, x 1 q = Vq cos <Pq 
V~Vg = psin6», 2/g=Zgsin0g 
dxdy — pdpdO. (B5) 

Then U 4 in Eq. (|B2|) becomes 

OO 27T 

I : = I p dp [ d0' e ^{VG) 2 +P 2 +2pi%co S {e'+<f, i:i -<t>%)) 



o o 



4^ 



U2e+p 2 + r -f ) 2 -p 2 4cos 2 (^)j 



(B6) 



where we moved to a new angle variable 9' — 9 — (f>ij. 
Using the expansion 

OO 

e -Xcos(e'+ a ) _ (-l) n I n (X)e ina e in8 ' (B7) 



71 — — OO 



in Eq. (|B6[) and keeping the n = term, the integration 
over the angle variable leads to: 

dff ^ 



--, (B8) 



where we used the formula J d9 
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c+b cos 6 ^/(c-b)(c+b) ' 

We shall proceed to p integration. Here we should 
note: In the case £ <C 1, the Lorentzian-like func- 
tion on r.h.s. of Ea. flBgp is sharply peaked around 
p and is well approximated for p > by Gaussian, 

,.h exp 



(p- 



4£2 



Therefore, the p integration 



leads to: 



u. 



o 




exp 



(P- 



4£ 2 



dp 



G ri i 



(l+£ 2 /A)(l + 8£ 2 /A)i/2 r 



x exp 



2% 



exp 




where the integration was carried out by the saddle-point 
approximation which is justified in the case £ 1. Since 
each vortex dynamics occurs in the range a A 13 <C 1 
because of the notion above Eq. fTo) in Sec|TT] , we may 

approximate Iq ( ^(f+aP/A) ) ~ ^ an< ^ ne gl ec t the contri- 
butions from I n (x) with n = 1,2. • • • in the region 
Then we reach 

4 



04 



^V2A 2 r 2 - 



*j\ — 



ff2 
7T 2 A 2 



iV 

-Vo 



U, 



exp 



2A 



(BIO) 



with use of (-^) ~ 1- At first the magnitude of the 
inter-vortices interaction looks very small (i.e., 0(£ 4 )), 
but, after scaling to make unity the inertial mass of each 
vortex, it becomes 0(92) (see Eq. ljP?)) ). It is interest- 
ing that the inter-vortices interaction energy in a multi- 
component BEC with no trap estimated in a different 
approximation (i.e., under the Abrikosov ansatz) is also 
proportional to -3- in the asymptotic region 31 1. 

Finally we note the remaining contributions, Uq + U2, 
in Eq. (|B2p . Their integration is quite simple and gives 
rise to 



u 2 =. 

92 

ttA 



92 
7T 2 A 2 

N 



2^(7- 



u 2 



2$ 2 



(Bll) 



10 



In the final expression, the first term is constant, giv- 
ing no contribution to the dynamics in Eq. (ll3[) . and the 
second one is a sum of single-particle contributions of 



0(t±) which renormalizes the 7-th line in Eq.®, giving 
no substantial contribution to Eq. (|13l) . 
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